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Abstract 



Following Feynman's successful treatment of the polaron problem we apply 
the same variational principle to quenched QED in the worldline formulation. 
New features arise from the description of fermions by Grassmann trajec- 
tories, the supersymmetry between bosonic and fermionic variables and the 
much more singular structure of a renormalizable gauge theory like QED in 
3 + 1 dimensions. We take as trial action a general retarded quadratic ac- 
tion both for the bosonic and fermionic degrees of freedom and derive the 
variational equations for the corresponding retardation functions. We find a 
simple analytic, non-perturbative, solution for the anomalous mass dimension 
7 m (a) in the MS scheme. For small couplings we compare our result with re- 
cent four-loop perturbative calculations while at large couplings we find that 
7 m (a) becomes proportional to ^fa. The anomalous mass dimension shows no 
obvious sign of the chiral symmetry breaking observed in calculations based 
on the use of Dyson-Schwinger equations, however we find that a perturba- 
tive expansion of 7m (a) diverges for a > 0.7934. Finally, we investigate the 
behaviour of j m (a) at large orders in perturbation theory. 
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I. INTRODUCTION 



Variational methods are widely used in many areas of physics but are not very promi- 
nent in field theory []J. This is due to the infinite number of degrees of freedom and the 
singular short-distance behaviour of relativistic field theories. A very successful application 
of variational methods in a non-relativistic field theory is provided by Feynman's treatment 
of the polaron @: after integrating out the phonon degrees of freedom and approximating 
variationally the remaining effective action by a retarded quadratic trial action one obtains 
the best approximation scheme which works for both small and large coupling constants. 
Detailed numerical investigations || have shown that Feynman's approximate solution devi- 
ates at most 2.2% from the true ground state energy for all coupling constants. It is therefore 
very attractive to apply similar techniques to problems in relativistic quantum field theory 
where there is much need for non-perturbative methods. In previous publications we have 
done that in the context of a scalar, super-renormalizable model theory ||. 

In this paper we present the first results obtained by applying polaron variational meth- 
ods to a realistic theory, namely Quantum Electrodynamics (QED) in the quenched approx- 
imation where electron-positron loops are neglected. While the actual coupling constant 
between electrons and photons a = e 2 /(Ait) ~ 1/137 is small enough to apply perturba- 
tion theory in most cases, there is enough interest to study the theory at larger coupling: 
first, the strong coupling behaviour of any physical theory is of interest in itself, second, the 
possibility of chiral symmetry breaking || demands an investigation at large a and, finally, 
bound state problems are inherently non-perturbative and involve powers of lnl/a ~ 4.92 
in radiative corrections. 

The extension of our methods to QED requires a formalism to include fermions and 
a treatment of the more severe singularities encountered in a renormalizable field theory 
rather than a super-renormalizable or non-relativistic one. We do this within the worldline 
technique which has recently experienced a revival ||. In this formulation, the degrees of 
freedom describing the electron are its bosonic worldline x^t), which is the four-dimensional 
analogue to the polaron trajectory, as well as a Grassmannian path (^(t) needed to describe 
the electron's spin |7j . Here t is the proper time which parametrizes the paths and runs from 
to T. The dynamics of the electron in an external vector field A^(x) with field strength 
F^x) are then described by the following worldline Lagrangian 

L = _^2 + iC .£ + L ±<x - e± . A {x)--F^{x)t% v . (1) 

Here k,q is an arbitrary parameter which may be used to reparametrize the proper time 
without changing the physics and \ is a Grassmannian (super-)partner of the proper time 
T. Note that the above action exhibits a well-known supersymmetry between bosonic and 
fermionic degrees of freedom ||. For further details about the application of the worldline 
formalism to QED we refer the reader to Ref. |§. 

The photon field may be integrated out exactly in complete analogy to the phonons 
in the polaron case, resulting in an effective action for the electron only 
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and G^ik) the gauge-fixed photon propagator. As described in Ref. ||, the electron prop- 
agator is obtained by carrying out a path integral over the degrees of freedom x^t) and 
Cn{t), as well as a weighted integral over the proper times T and \i an d finally identifying 
the Grassmannian variable V = £(0) + C(T) with the Dirac matrix 7: 
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Here M is the bare mass and T>x contains an integration over the endpoint x = x(T). 
Note that we have divided and multiplied by the path integral for the free theory, so the 
bare propagator may be obtained by just ignoring the last line. For non-zero couplings, of 
course, the path integrals in the last line cannot be performed; these we shall approximate 
variationally in the next section. 



II. VARIATIONAL APPROACH 

Feynman's variational principle has its root in Jensen's inequality for convex functions 
applied to exp(— Se), where is a Euclidean action. In Minkowski space and/or for 
complex actions the variational principle remains valid, however it becomes a stationary 
principle rather than a minimum principle. To be more precise, the path integral over 
bosonic and fermionic paths obeys 
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where <...>( indicates an average involving the weight function e tSt in the relevant func- 
tional integral and S t is a suitable trial action. Note that corrections to this variational 
approximation may be calculated in a systematic way and that, furthermore, to first order 
in the interaction (i.e. to order a) the relation is in fact an equality if St reduces to the free 
action for small couplings. 

For the trial action required in Eq. (|^) we choose a general retarded quadratic action 
which is a two-time modification of the free action in Eq. (Q) 
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Here the variational parameters are contained in the retardation functions gi(cr) for bosonic, 
fermionic and spin-orbit interactions; these are even functions of a = t\ — £2 and they 
become identical for a supersymmetric trial action []. The variational principle 'adjusts' 
these functions in order to compensate for the fact that the true effective action ([|) is not 
quadratic in the variables x(t), ((t). Feynman's polaron result was obtained by taking a 
specific Ansatz for the retardation functions but here we leave their functional form free. 
This is because one expects that the correct short-time behaviour of these functions is 
much more important for a renormalizable theory like QED than for the polaron problem 
which does not exhibit any ultraviolet divergences. Indeed one finds that for small a the 
"best" ^^(cr) behaves like y/a, In a and 1/cr in the polaron, super-renormalizable and QED 
case, respectively. The "tilde" over the trial action indicates that it includes the boundary 
terms already present in Eq. (f|) and that we are using "momentum averaging" [ffl. These 
terms, involving the external momentum p and the Grassmann variable T, are multiplied by 
additional variational parameters Ai and A2, respectively. They provide additional freedom 
to modify the strength of the boundary terms. We have allowed this freedom because of our 
experience in scalar relativistic field theory J3J, where the variational parameter Ai turned 
out to be essential for describing the instability of the Wick-Cutkosky model. 

Since the trial action @ is at most quadratic in x(t) and ((t) it is possible to evaluate 
the various averages required in Eq. (|^) analytically. A particular simplification occurs if 
one restricts oneself to p 2 = M 2 , where M is the physical (i.e. pole) mass: as discussed 
in [|J, the divergence of the propagator on its mass shell results from a divergence of the 
integral over the proper time T. Indeed, the variational approximation fl5|) results in an 
electron propagator [see Eq. (^)] which has the form 

GT{p) = e^w jT dT J d X exp [-M 2 + p 2 (2X - A 2 )] } (7) 

.exp\-—(n[A B ]-n[A F )+V\p&,iJ? F )) + F( X ,r;T;p)j , 
I v ' J r=o 

where A (which is defined below), the f2's and V are T-independent and the function 
F(x, T; T; p) is subleading in T. The latter therefore contains information relevant for the 
wavefunction renormalization of G^ip), and not the pole structure. We leave the discus- 
sion of this function, which also contains the entire X an d T dependence (and hence the spin 
structure of the propagator), for a future publication as it is not required for our present 
investigation. 

From Eq. (|7]) we see that the bare and physical mass are related through 

M 2 =M 2 (2A-A 2 )-2( n[A B ]-n[A F ] + V\jj* B ,t4] ) . (8) 



1 We have explicitly separated out the free action in Eq. (g), which could have alternatively been 
added into the second term by adding 5(a) / (2ikq) to each of the gi(a)'s. This way our retardation 
functions contain no distributions. Also, in the supersymmetric limit our trial action could be 
written in the explicitly supersymmetric notation of Ref. as Ikq J dtidt2 J d9\d02 g{T\2) DX\ ■ 
DX2, with a single retardation function. 
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We have labeled this relationship Mano's equation as K. Mano first applied polaron tech- 
niques to a scalar relativistic field theory []T0| | . Note that, on mass shell, the variational 
equations resulting from Eq. (|5|) are equivalent to demanding stationarity of Mano's equa- 
tion. 

The nomenclature in Mano's equation corresponds to that introduced in Ref. [|J: f2L4#] 
and fipAp] originate from contributions (bosonic and fermionic, respectively) of the terms 
in Eq. (|5|) involving Sq and St only. They are the analogue to the kinetic term in variational 
quantum mechanical calculations, while the analogue of the contribution from a potential 
term (explicitly proportional to the strength of the coupling) resides in V. 

Similarly to Ref. [[|, it is useful to express the retardation functions in terms of the 
variational "profile functions" Ai(E) and the "pseudotimes" [if (a) ,i = B,F defined by 
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respectively Q. Furthermore, it is convenient to define A = X\/Ab(0). Indeed it turns out 
that the averages in Eq. (|5|) can be directly expressed in terms of these quantities. The 
kinetic terms in Eqs. (0) and (|8|) become 
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where d is the spacetime dimension d = 4 — 2e. This is identical to the result in Ref. |4j] if 
d = 4 and Kq = i (i.e. the Euclidean formulation) are taken. The specific properties of QED 
are encoded in the "interaction" term V which, with V = V\ + V2, reads 
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Note that by fi 2 {a) we mean -^fi 2 (a) (and not -^fx(a) ), the function E(k,a) is defined 
to be E(k,a) = exp {i \k 2 ji 2 B {a) — 2Xk - pa] /(2k )} and of course p 2 = M 2 . The fermionic 
contributions, both in the 'kinetic term' Qj? as well as in Vi, appear with an opposite sign 
to the bosonic contributions. The reason for the separation of V into two pieces will become 
apparent below. 



2 It is easy to show from Eq. ( |j~0| ) that, as long as Ai(0) and Aj(oo) exist, for both asymptotically 
large and infinitesimally small times a the functions fi 2 (a) become proportional to a; hence the 
label "pseudotime" . In the free case one has M E ) = b/^iV) = a,X = = 0. The other 
variational function Ago(E) and parameter A2 are linked to the spin structure of the propagator 
and therefore do not show up in Mano's equation. 
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By construction Mano's equation is stationary under variation of the parameters. It 
is important to note that we have not demanded the various retardation functions g B ,F 
(as well as gso-i which only plays a role for the residue) to be identical (before variation). 
Had we done so, the resulting profile functions A B and Ap would have also been identical, 
the pseudotimes fi 2 B F would have been one and the same and hence Q[.Ab] — fiLAp] as 
well as V\ would have vanished. The absence of a 'kinetic' contribution would have been 
fatal to the variational principle as this contribution provides the restoring 'force' to the 
potential V. On the other hand, closer examination of V\ reveals that (x 2 B ^ ftp is also 
dangerous: The contribution of each of these terms is quadratically (UV) divergent if the 
dimensional regularization is replaced by a momentum cutoff. This may be checked by either 
directly substituting the small a limit of /xf (cr) into V\ or by noting that for scalar QED, 
where the Grassmannian path integrals are absent, the remaining contribution from jx B gives 
rise to the quadratically divergent one-loop diagram of that theory. It is the combination 
(fi 2 B ) 2 — (ftp) 2 which displays the usual logarithmic UV divergence of QED. Although at 
leading order in the coupling we are guaranteed to reproduce the correct perturbative result 
[see Eq. (|5|)], at higher orders the cancellation of these quadratic divergences is ensured by 
the super symmetry. To summarize, on the one hand the trial action cannot be restricted 
to contain only supersymmetric terms but on the other hand allowing non-supersymmetric 
terms may destroy the renormalizability of the theory. 

The way out of this predicament is provided by the variational principle itself: although 
it is unavoidable that the trial action breaks supersymmetry, the actual solutions to the 
variational equations may in fact be nearly supersymmetric. That this indeed turns out to 
be the case may be seen by recognizing that V\ is the most singular part of the interaction 
whereas V2, which involves only bosonic contributions and is the only source of supersymme- 
try breaking, is similar in structure to the scalar super-renormalizable model studied before 
||1 1|| . Divergent contributions in the limit e —>■ to the variational equations are solely de- 
termined by V\. Therefore, the divergent contributions to A B (E) and Ap(E), and hence to 
ji B and fip, are identical. 

In this paper we confine ourselves to studying this divergent structure and so it is suf- 
ficient to set A B (E) = Ap(E) = A(E). The corresponding variational equation becomes, 
after performing the fc-integration in Eq. (fT~2|) . 
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Note that here we have now also dropped the subscript on the pseudotime as it is no longer 
relevant and we have defined 

a /2Tci\ e 3 - 2e e ^ 3a 
° e ~ ttUo J (l-e)(2-e) "~ ' 2vr ' [Lb) 



Since fi 2 (a) — > a for small a one sees that the u-integral in Eq. flT4" ) would diverge for e = 0; 
this just reflects the 1/a behaviour of the retardation function in Eq. (0) as was discussed 
before. The crucial difference between super-renormalizable and renormalizable theories 
therefore is that for the latter ones the variational equations themselves are UV-divergent. 
In this way the divergent structure of higher-order diagrams is effectively summed up. 
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We may now simplify V by making use of the above "asymptotic" supersymmetry. The 
only remaining contribution is that of V 2 = X 2 M 2 W2 which becomes, after carrying out the 
integration over the momentum k, 



\--e)(l-e) 2e f°° da r 1 du . M , . ft X 2 M 2 a 2 » 
10 lfi 2 (a)} jo u \ 2k n 2 (a) 



where the u integration arises from an exponentiation of the photon propagator in Eq. ([13]) 
in a similar way as in Ref. [|J. With this, the variational equation for A in this asymptotic 
limit becomes 

A = 1 - ^ 2 W 2 ) . (17) 



III. MASS RENORMALIZATION 



Renormalizability of (quenched) QED means that all divergences can be collected in 
the mass and wave function renormalization constants. In the present investigation we 
concentrate on the mass renormalization constant in the MS scheme, Zj^ s , defined via 
M = M v where M u is an intermediate mass scale. In this scheme it has the perturbative 
expansion 



7 MS 
J M 



b u a 



+ 



b 2 2 


61a 





_e 2 " 


e 





+ 



(18) 



where it is known from perturbation theory [|T!| that the expansion coefficients bij are pure, 
i.e. mass independent, numbers. Furthermore, the renormalization group provides relations 
between many of these coefficients; at order n in perturbation theory only the coefficient b\ n 
contains new information. This is encapsulated in the solution of the renormalization group 



equation for Z^ s , namely 
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where j m (a) is the anomalous mass dimension of the electron |T3|| . In perturbation theory, 
7 m (a) can be extracted from perturbative QCD calculations, which have been performed 
up to 4-loop order. One obtains 70 = 3/2, 71 = 3/16 |L4j], 72 = ^ = 2.0156 and 73 = 



1 

128 



^ + 336C(3) 



As the variational 



= -4.3868 U). 

calculation is applicable for arbitrary values of the coupling, compar- 
ison to perturbation theory provides a useful guide to its utility. As mentioned before, to 
first order in the coupling the calculation is guaranteed to be exact as long as one has used 
a trial action which can reduce to the free action in the limit a — > 0. A genuine test of the 
variational scheme is only obtained by comparing the coefficients in higher order. It should 
be noted that this test is much more demanding than in the polaron case where one can only 
compare the numerical value of the second-order coefficient for the energy: here, in addition, 
one tests the e-dependence of this coefficient and also whether it is mass-independent as it 
should be in the exact theory. 
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In order to know and V at second order in a one requires the variational parameters up 
to first order in a. These may be obtained by inserting the zeroth order results /i 2 (o") = o 
and A = 1 into the variational equations for the profile function ( |i4|) and A (|T7D. The 
solutions then need to be substituted back into Q [Eq. ([TT|) 1 and V 2 [Eq. fllBD l. Having 
done this, Z^ s may then be extracted from Mano's equation, yielding = 9/32, which 
is correct, and fe^f = 0, which should be compared to the exact value of b\2 = —3/64. As 
in the Wick-Cutkosky model, the A-variation is of crucial importance: for example, fixing 
A = 1 would give a wrong result for 6 2 2 and a logarithmic mass-dependence for 6 12 . 

It is possible to develop the perturbative expansion of the variational result further, with 
the result that no mass dependence in the coefficients appears even at higher order. Indeed, 
it turns out that it is in fact possible to obtain the full analytic expression for the anomalous 
mass dimension in the worldline variational approximation. We shall sketch the derivation 
below, leaving the technical details for the Appendix to this paper. 

To begin with, we first drop the mass term in the variational equation (|T4]) since it only 
affects long-distance physics and not the ultra-violet behaviour contained in Z^ s . Then 
we change variables from a to y = c e (u 2 ay, and equivalently for E. This has the effect of 
making the system of integral equations (|H],|14],|]1|) independent of the coupling. We write 
these explicitly in the Appendix, where it is shown that if, for small e, the pseudotime has 
the form 



a*V) [ My) . nf 

= exp h (J(( 



(20) 



a 

then the anomalous mass dimension may be written in terms of this function Uo(y), i.e. 

lm = T^R ' (21) 

where v(y) = y Lo' (y) and y is determined by the implicit equation 

Vo = ^ e Mm) • (22) 

On the other hand, it can also be shown (see the Appendix) that the variational equation 
for the pseudotime translates into an equation for u> (y), i.e. 

nl 1 -v{y)\cot 



, . , . (23) 

y 2 L L2 -'J 1 ; 

This equation must in general be solved numerically. However, we note that the calculation 
of the anomalous mass dimension in Eq. fl2T|) only requires knowledge of the function v(y) 
at y = yo. Furthermore, it is remarkable that, at this value of y, the combination e^^/y 
is precisely the combination that is fixed in terms of the coupling constant [see Eq. (p2])j. 
Hence, at y = yo, the the L.H.S of Eq. (|23|) may be written in terms of a while on the R.H.S. 
we can eliminate v(yo) completely in terms of 7^ r by making use of Eq. fl2"I|). One is left 
with a simple implicit algebraic equation for the anomalous dimension 

l a = (l + ^tanf 7 / 27 ^ , (24) 

without ever having actually solved the variational equations themselves. Eq. ([24]) is the 
main result of this paper. 
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IV. DISCUSSION 



When expanded in powers of a, Eq. fl2~4] ) immediately yields 

var/ \ 3a 9 2 /a\ 3 27 2 /a\ 4 243 2 / vr 2 \ /a\ 5 ^. 6 . 
7 " (0) = 2 7 - 32^ U + 32* ) - 128^ I 1 - 20 j UJ + ° (a » ■ (25) 

which may be compared to perturbation theory. Numerically the values of the coefficients 
are different but of the same order of magnitude as the exact perturbative results. Note, 
however, that this comparison is not particularly meaningful: the variational result is an 
approximation which is valid at all a. It need not have the same, or even approximately the 
same, perturbative expansion in a as the exact result. It should, however, be numerically 
similar. In Fig. [I] we plot the variational result as a function of the coupling and compare 
it to perturbation theory up to 4-loop order. For a £ 1 the 3- and 4-loop anomalous 
dimensions start to deviate so much from each other that one cannot trust either of them. 
Also shown is the result up to 5 loops, where the 5-loop coefficient has been estimated from 
Pade approximations to the perturbation theory (see Eq. (2.12) of Ref. ||16|| , which needs 
to be adapted to QED with n/ = flavours; one finds <-y| >ade = 3.848 ). Clearly this 
does not significantly extend the numerical validity of the perturbative result. In short, the 
variational estimate for 7 m is roughly in agreement with (albeit apparently a little below) 
the perturbative result in the region where the perturbative result can be trusted. 




FIG. 1. Anomalous mass dimension 7 m as function of the coupling constant a in quenched 
QED. The variational result (^4|) is shown as a solid curve while the solution from the 
Dyson- Schwinger equations in rainbow approximation is indicated as a dot-dashed curve. The 
curves labeled "n-loop" show the result up to n-loop perturbation theory. Finally, the Pade esti- 
mation of the 5-loop result is also shown. 
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Also shown in Fig. [I] is the only other easily available non-perturbative result for 7^ s 
based on the use of dimensionally regularized Dyson-Schwinger (DS) equations in "rainbow 
approximation" within the Landau gauge. This may be obtained by adapting the discussion 
in Ref. |17j to finite Mq, with the result that 7^ s = 1 — Jl — 3a/n (which is the same as 



derived by Miransky [y| using a hard momentum cutoff). We see that this result deviates 
from perturbation theory in a region where, at least numerically, perturbation theory still 
appears to converge. Above a = 71/ 3 = 1.047 the DS result becomes complex, this value of 
the coupling constant coinciding with the coupling a cr at which the onset of chiral symmetry 
breaking takes place in those calculations. This is in contrast to the variational result which 
remains real for all values of the coupling and in fact has the strong coupling limit 

7 r(«) a -^\^-\ + °{-^=) ■ ( 26 ) 

Further investigations are necessary to clarify the absence of any obvious sign of chiral 
symmetry breaking in the variational result for 7^ s (a) 0. Indeed, in order to investigate the 
issue of dynamical chiral symmetry breaking, it would seem to be more straightforward, at 
least conceptually, to set M on the right hand side of Mano's equation (||) to zero and to 
see if the variational equations can be satisfied in this case (for a finite physical mass M). 
This, however, goes considerably beyond the scope of this paper: we have merely calculated 
Zf\ s = M /M u (or, more precisely, 7™ S ), which meant that we could simplify the calculation 
by i) restricting ourselves to considering the supersymmetric massless limit of the variational 
equations in Section II and it) only taking into account the most divergent contributions (as 
e — > 0) to the variational equations, as well as to W%, in Section III and the Appendix of this 
paper. In a full calculation of the R.H.S. of Mano's equation (i.e. the additional calculation 
of the finite renormalization M v /M) these two simplifications should not be made 0. In 
other words, even if Z^ s 7^ 0, dynamical chiral symmetry breaking can still occur if M v /M 
vanishes for finite M. 

It is interesting to note, however, that there are also some strong similarities in the 
analytic structure of the variational and DS result. A perturbative inversion of Eq. (|24]) , 
i.e. the expansion 7™ r (a) = J2™=i c n a n , has a finite radius of convergence due to a branch 
cut in the complex a plane. The position of this cut, and hence the radius of convergence, 
can be determined most easily by searching for the value of a at which Eq. (E3) has two 
solutions for 7™ which are infinitesimally close to each other. This amounts to demanding 
that Eq. (J33j) is satisfied and at the same time the derivative of its R.H.S. vanishes, i.e. 

= cot(^) + -^/sin>(-^U . (27) 



3 The reader should note that the issue of dynamical chiral symmetry breaking in a dimensionally 



regulated theory is a notoriously subtle problem; see Ref. 17]. In particular, it was shown there 
that if four dimensional quenched QED breaks chiral symmetry above a critical coupling than the 
dimensionally regularized theory will break it for all couplings at finite e. 

4 This situation is analogous to what is the case in perturbative calculations, where anomalous 
dimensions of operators are far easier to calculate than finite contributions. 
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One finds that a con = 0.7934, which is not too different from the radius of convergence of 
the DS result It is not clear whether this similarity between a cr and a con is accidental 



or not. 

In connection with this, it is interesting to note that for large n the behaviour of the 
expansion coefficients c n in both the variational result as well as the DS result are rather 
similar: 



cx ,-, /ri sin 



5tt\ 3tv 
a + — n 



7 J 7 



(28) 



where numerically (3 ~ 1.38, a ~ 2.3 x 10 -3 and b « —8.27 x 10~ 2 . For the DS result one 
obtains j3 = log(2y / 7r) = 1.27 and the sine function is absent. It is the sine function in 
the variational result which is responsible for placing the branchpoint (which, for the DS 
result, is on the positive real axis) into the complex plane. Furthermore, it is remarkable 
that the large-a limit of |7 m (a)| obtained in Eq. ( p6|) is almost the same as for the DS result: 
\j™(a)\ - 1.0 vs. | 7 ° s («)| - 0.98^. 

It should be pointed out that a finite radius of convergence of the perturbation expansion 
is not what one generally expects from calculations of large orders of perturbation theory 
using the methods of Lipatov and others p0| . Rather, the factorial growth of the number 



of diagrams at ia th order in perturbation theory tends to lead to a vanishing radius of 
convergence. As has been observed elsewhere [|J, it can be shown that the variational 
calculation contains (pieces of) all possible Feynman diagrams at any order in perturbation 
theory. One concludes, therefore, that at n th order in perturbation theory there are either 
strong cancellations between diagrams in the variational calculation or that 0(n\) of them 
give a vanishing contribution. 



V. SUMMARY AND OUTLOOK 

We have applied polaron variational techniques to quenched QED in 3 + 1 dimensions 
and obtained, within the MS scheme, a remarkably simple expression for the anomalous 
mass dimension valid for arbitrary couplings. The approach has considerable advantages 
over other techniques in that it automatically maintains gauge invariance, as well as the 
requirements of the renormalization group, and corrections can be systematically calculated 
(as has been done in the polaron case pl|). Furthermore, we have shown that the numerical 
results for / ~f m are rather reasonable at small coupling and that at large couplings the per- 
turbative expansion of this quantity fails in a way similar to rainbow DS results. It would 
be interesting to compare to DS calculations which go beyond the ladder approximation, 
thus decreasing the strong gauge dependence inherent in that approximation. Furthermore, 
variational calculations with more general trial actions could give an indication whether this 
analytic structure is robust, thus indicating possible large cancellations between diagrams at 
high order in the perturbation theory of quenched QED, or whether this structure is just an 
artifact of the particular trial action used in this paper. Finally, we note that the calculation 
of physical observables or application to bound state problems also seem feasible within the 
variational worldline approach developed here. 
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APPENDIX 



In this Appendix we provide some of the technical details which enter into the derivation 
of the variational approximation to the anomalous dimension. To begin with, we shall scale 
the trivial a dependence out of /x 2 (cx) and define the reduced pseudotime s(a) as 



a s(a\ 



(Al) 



As argued in the main text, mass terms can be dropped for the calculation of the mass 
anomalous dimension. A perturbative evaluation of the variational equations (|14]) and (|10D 
for M = then shows that profile function and reduced pseudotime have an expansion in 
powers of and a e , respectively: 



(2 \ ne 
^ ,*{*) = i + 5> n («/V) n 
^ J n=l 



(A.2) 



One finds 

A,: 

S2 : 



c e T(e) COS 



1 (Ce 

2 U 

1 (Ce 

2 Ve 



1 



1 + e 

2 1 
l + 2e 



3a 1 
27r e 



T(l + 2e) cosfevr) 



1 + 



e->o 1 /3al \ 2 

* 2 V 2tt e J 
r 2 (l + e) 1-e 



3a 1 
"2tt e 



(A.3) 



~^ 2 



cos(eTr) J T(l + 2e) 1 + e 
This suggests that perhaps the leading e-behaviour of the coefficients is 

. -o (-1) 



3a 1 
2^7 



.4, 



e->o J_ / 3a I 
' n\ V 2tt e 



nl 



3a 1 
2vr 7 



(A.4) 



In view of these results we rewrite all equations in terms of the dimensionless quantities 

,2\ e 



V 



c f \vo 



z 



E 



(A.5) 



where we have also rescaled by c e (which is linear in the coupling) because the mass scale v 
always appears in the combination a v 2e in dimensional regularization. In a similar way, the 
variational parameter A almost always appears in the combination AM, so it is convenient 
to define the dimensionless combination 



' 2k u 2 " 
i\ 2 M 2 



(A.6) 
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With these definitions, the massless variational equation (|14]) for A(E) may be brought 
into the form 

a , \ -, 1 f°° , cosiy/z) 1 ^ . t . 

A(z) = 1 + - / dy ) y 1 ' . A.7 

The reduced pseudotime [ i.e. the rescaled version of Eq. (|ToD ] is now given by 

2 1 f 00 , I f z\ 1/e l-cosiy/z) 1 /" /A s 

s(y) = — dz-(-) ) y ! } , A.8 

7r e jo z \y J A{z) 



while the rescaled potential W2 becomes a function of a e alone: 

~(y/ae 



W 2 (a £ ) = - / } . ^ / — e + (1 - e) u] exp 

and hence the variational equation ( |TTD for A becomes 

1 



i/ e JL 



s(y) 



(A.9) 



1 + 2W 2 {a e ) - 2ea e W! 2 {a t ) . (A.10) 
A 



The anomalous mass dimension in the MS scheme may be defined [ see Eq. (|19|) ] through 

a d 

J A1 



lm = -H'zbfo 2 " • (A - n) 



Note that this equation is correct independently of whether Z M has been calculated in the 
MS scheme or whether it is defined through Mano's equation by Zm = M /M. It is because 
of this fact that we can derive a nonperturbative expression for 7 m in the MS scheme, even 
though this scheme is usually only used within the context of perturbation theory. 

As the a dependence of Zm now only enters through the variable a e (and of course 
implicitly through the variational parameters), it is not surprising that we may use the 
variational equation ( A.10|) for A to simplify j m . Indeed, differentiating Mano's equation 
with respect to the coupling gives 

The first term is zero because of the variational equation for A, da e /da is just a e /a and by 
substituting the variational equation for A into Mano's equation we find 

Z 2 M = A[l-2eAa e W 2 '(a e )]. (A.13) 

Hence the anomalous mass dimension is just given by 

7m = , im . . (A.14) 

In order to proceed further, we need to evaluate W 2 (a e ). In general one would need to 
do this numerically, however fortunately in Eq. ( |A.14j) only the small-e limit is required. Let 
us assume that the reduced pseudotime may be written as 
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s(y) = exp 



w(l/,e) 



(A.15) 



where u)o(y) = lim 6 _, u(y,e) is finite. This is supported by the perturbative results given 
in Eqs. (|A.2| ) and (|A.4| ) and we shall show that this holds in general when we solve the 
variational equations below. In this case the exponential in Eq. (|A.9 ) has the argument 



(A.16) 



If the term in brackets is larger than one, this argument will become arbitrarily large (and 
negative) in the limit e — > 0, hence it will lead to a vanishing contribution to the integral. If 
the term in brackets is smaller than one, however, the argument goes to zero, the exponential 
factor in Eq. ( |A.9|) may be replaced by unity and the integral over u may be performed, 
yielding 



W 2 (a e) 



1 fVO 



2e 



dy exp 



u(y,e) (1 



- € 



where y is given by the equation 



y e Mvo) = lima 



(A.17) 



(A.18) 



We have assumed here that ye^ ^ is an increasing function of y, which will turn out to be 
the case. The leading term in Eq. ( |A.17|) may be obtained by integration by parts, with the 
result 



W 2 (a e) 



2co'(y , e 



exp 



[l-e] 



(A.19) 



We also require the derivative of this function, which is most easily obtained by direct 
differentiation of Eq. ( |A.17|) : 



Wi(a e ) ^ lexp 



w(j/o,e) 



'l-e) 



-^(yo,e) 



1 + y u'(y ,e) 



(A.20) 



Substitution into the variational equation for A yields 

1 + y u'(y ,e) 



A 



£^0 



w'(yo, e) 



and hence 



2e\a f W'(a f 



exp 



so that the anomalous dimension becomes 

Vou'oiyo 



Vo u'(y ,e) 



7r7 



i - you' (yo) 



(A.21) 



(A.22) 



(A.23) 
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We stress that only the last line is exact while the previous ones have correction terms for 
finite e. In particular, the calculation of Z M = M /M (as opposed to Z^ s = M /M u ) would 
require these additional terms and hence the result, unfortunately, does not shed light on 
whether Zm could in fact be zero for finite M, which would signal chiral symmetry breaking. 
Note that the limit e — > in Eq. ( |A.18| ) needs some care: naively, one would conclude from 
the definition ( A.6 ) that the R.H.S. equals lim e ^oc e = 3a/(2n) but Eq. ( |A.21| ) shows 
that the variational parameter A vanishes like exp(— Uo(yo)/^) arid therefore also gives a 
contribution: 



2/o e 



3a 
2^ 



,2u} (y ) 



(A.24) 



this being the result (0) quoted in the main text. 

It now remains to calculate the function u>o(y). The arguments used to derive the ap- 
proximate expression for W^&e) m Eq. ( |A.17|) are more difficult to apply to the variational 
equation (|A.7|) for A(z) and the definition ( |A.8|) of s(y) because of the rapidly oscillating 
trigonometric functions appearing in their integrands. We shall therefore adopt a more sys- 
tematic approach at this stage and note that it is possible to write these equations in a 
differential form. Consider, for example, an integral of the type 



dyf(y) cos 



r 



(A.25) 



By changing integration variable to y 1 ^ and Taylor expanding the function f(y) we may 
carry out the integration term by term by making use of the integral 



dy y q 1 cos y 



T(q) cos ( q- 



71 



Hence we obtain 



h{z) = £/ (n) (o) 



yll+l 



n=0 



n + 1 



r [1 + (n + l)e] cos 



(n + l)e- 



71 



(A.26) 



(A.27) 



dyf(y) =: lc(eD z ) / dy f(y) 



This expression may be resummed, by defining the dilatation operator D z = z4-, into the 
compact form 

I t {z) = T{l + eD z ) cos(|eA 

Hence the variational equation for A(z) becomes 

A{z) = 1 + - lc {tD z ) [ Z dy 
e Jo 

and in a similar way we can rewrite Eq. ( |A.8j ) as 

S (V) = -ri^T 7~7TT "TTT • (A.30) 



(A.28) 



(A.29) 



1 + eD y i c (eD y ) A{y) 

Inverting Eq. ( |A.30|) and substituting into Eq. ( |A.29| ) eliminates the profile function A(z) 
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{l + eD y ) lc {eD y )s{y) 



1 [V 

1 + -Jc(eDy) dx 



s{x)\ 



l-e 



(A.31) 



Finally, one can eliminate the integral by operating with D y on both sides of this equation, 
so that 



D 



1 



y (l + eD y ) lc (eD y )s(y) 



-Jc{eD y 



y 



isiy)] 1 - 



(A.32) 



This equation may be solved systematically by defining s(y) in terms of the function u(y, e) 
[ see Eq. (|A.15[) ] and by making the Ansatz that u(y,e) has a power expansion in e 



^0,e) = uj {y) + eu x {y) + ... 



(A.33) 



The crucial observation is that repeated application of the dilatation operator on an expo- 
nential of the form of Eq. ( |A.15| ) results in (eD y ) n exp (uj/e) = \{yuj') n + 0(e)] exp (u>/e) so 
that at leading order in e, for any function F(eD y ) acting on exp(±u;/e), we have 



F(eD y ) exp(±co>/e) ^» F(±yuj' ) exp(±c<j/e) . 
Applying this relation to Eq. ( |A.32| ) provides the following equation for u>o(y) 

lc (yuj' )e-"° . 



(i -yu'o)ic(-yu'o) 



(A.34) 



(A.35) 



which is e-independent, justifying the Ansatz ( |A.33 ) a postiori. This equation may be 
simplified considerably by making use of the reflection formula T(z) T(l — z) = ir/ sin7rz for 
T-functions. By defining v (y) = y uj' (y), we then find 



71 



1 — v (y)] cot 



7T 



v(y) 



(A.36) 



which is Eq. (^3|) in the main text. Together with the boundary condition uq(0) = [ i.e. 
/x 2 (<t) — * cr for a — > 0, as discussed in the footnote below Eq. fliTf ) ] the first-order nonlinear 
differential equation (|A.36|) determines the function u)o(y). Remarkably, as shown in the 
main text, it is not actually necessary to solve it in order to obtain the anomalous mass 
dimension j m . It is also interesting to note that due to the reflection formula all T-functions 
have disappeared, which has the consequence that in a perturbative expansion of 7™ r no 
Riemann (^-functions, but only powers of tt, occur. 
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